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Abstract 


The ground beetle genus Agonum Bonelli, 1810 is a large genus of the tribe Platynini with many species that show high amounts of 
intraspecific variations, making a correct identification challenging. As part of the German Barcode of Life initiative, this publica- 
tion provides a comprehensive DNA barcode library for species of Agonum that are reported for Germany. In total, DNA barcodes 
from 258 beetles and 23 species were analysed using the Barcode of Life Data System (BOLD) workbench, including sequences 
from former studies and 68 newly-generated sequences. The neighbour-joining analyses, based on K2P distances, revealed distinct 
clustering for all studied species, with unique Barcode Index Numbers (BINs) for 15 species (65%). BIN sharing but distinct clus- 
tering was found for three species pairs: Agonum micans/Agonum scitulum, Agonum impressum/Agonum sexpunctatum and Agonum 
dufischmidi/Agonum emarginatum. The given dataset and its analysis represent another important step in generating a comprehensive 


DNA barcode library for the ground beetles of Germany and Central Europe in terms of modern biodiversity research. 
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Introduction 


Species identification represents a pivotal component for 
biodiversity studies and conservation planning, but rep- 
resents a challenge for many taxa when using morpholog- 
ical traits only (e.g. the correct identification of juveniles 
or larval stages). As a consequence of tremendous techno- 
logical advances in molecular biology during the last 20 
years, molecular data have become increasingly popular 
in species identification. In this context, DNA barcoding 
represents the central component in the modern diagnostic 
toolbox of molecular biodiversity assessment studies and 
taxonomic research (e.g. Hebert and Gregory 2005; Krees 
et al. 2015; Hajibabaei et al. 2016; Miller et al. 2016). For 
animals, a 658 base pair (bp) fragment of the mitochondrial 


cytochrome c oxidase subunit I (COI) gene has been select- 
ed as a DNA barcode (Hebert et al. 2003a, b). The utility of 
DNA barcoding relies on the assumption that genetic vari- 
ation within a species is much smaller than variation be- 
tween species. Barcode sequences are typically deposited 
in the international Barcode of Life Data Systems database 
(BOLD; http:/Awww.boldsystems.org). This public data- 
base acts as the central core data interface and repository 
that allows researchers to collect, to organise and to analyse 
DNA barcode data (Ratnasingham and Hebert 2007). In ad- 
dition to various analytical tools implemented in the BOLD 
workbench, DNA barcodes can be analysed using the Bar- 
code Index Number (BIN) system that clusters DNA bar- 
codes to produce operational taxonomic units that closely 
correspond to species (Ratnasingham and Hebert 2013). 
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Nevertheless, it should be noted that various effects 
can seriously limit the use of DNA barcoding or mito- 
chondrial markers in general, in terms of molecular spe- 
cies identification (e.g. Will and Rubinoff 2004; Will et 
al. 2005; Rubinoff et al. 2006; Collins and Cruickshank 
2013; Duran et al. 2019). For example, closely related 
but distinct species, as well as species that hybridise, 
may share identical haplotypes (e.g. Sota and Vogler 
2002; Takami and Suzuki 2005; Andujar et al. 2014). 
In contrast to this, complex phylogeographic histories 
(e.g. Schoville et al. 2012; Faille et al. 2015; Weng et al. 
2016) and incomplete lineage sorting effects (e.g. Sota et 
al. 2005; Zhang et al. 2005; Weng et al. 2020) can contort 
the mitochondrial variability of the studied organisms, 
generating para- and polyphyly in phylogenetic trees 
(Funk and Omland 2003; Ross 2014; but see Mutanen 
et al. 2016). Heteroplasmy (Boyce et al. 1989) or the 
presence of mitochondrial pseudogenes (Hazakani-Co- 
vo et al. 2010; Maddison 2012) can affect a successful 
amplification of the target fragment. Finally, maternal- 
ly-inherited endosymbionts, such as the alpha-proteo- 
bacteriae Wolbachia, may cause a linkage disequilibrium 
within mtDNA in arthropods and result in a homogenisa- 
tion of mtDNA haplotypes (e.g. Hurst and Jiggins 2005; 
Kolasa et al. 2018; Kajtoch et al. 2019). However, a vast 
number of studies across a broad range of taxa demon- 
strate the efficiency of DNA barcoding as the method of 
choice in molecular species identification (e.g. Raupach 
et al. 2014; Schmidt et al. 2015; Moriniére et al. 2017; 
Schmid-Egger et al. 2019). 

Not surprisingly, the build-up of comprehensive 
DNA barcode libraries represents a pivotal task (e.g. 
Brandon-Mong et al. 2015; Curry et al. 2018; Weigand 
et al. 2019). In Germany, the German Barcode of Life 
initiative (GBoL; www.bolgermany.de) aims to assess 
the genetic diversity of all animals, fungi and plants of 
Germany. Numerous studies provided first comprehen- 
sive DNA barcode libraries for various insect taxa (e.g. 
Hausmann et al. 2011; Hendrich et al. 2015; Hawlitschek 
et al. 2017; Havemann et al. 2018). In the case of ground 
beetles (Carabidae), a number of previous studies start- 
ed to build up a comprehensive DNA barcode library for 
these routinely-used biological indicators, including the 
genera Bembidion Latreille, 1802 (Raupach et al. 2011; 
Raupach et al. 2016), Amara Bonelli, 1810 (Raupach et 
al. 2018) and Notiophilus Dumeéril, 1806 (Raupach et al. 
2019). Detailed studies of many other important genera, 
however, are still missing. 

Within the Carabidae, the genus Agonum Bonelli, 
1810 is a species-rich taxon of shiny black or metal- 
lic ground beetles that imitate large representatives of 
Bembidion Latreille, 1802 (Lindroth 1986) (Fig. 1). Bee- 
tles of this genus are commonly found in forested and 
open habitats adjacent to freshwater throughout Europe, 
the eastern Palearctic area and North America (Liebherr 
and Schmidt 2004; Bousquet 2012). They are generalist 
predators that feed on a wide range of small arthropods, 
for example, springtails, aphids or midges (e.g. Griffiths 
et al. 1985; Bilde and Toft 1994; Hannam et al. 2008). 
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A first detailed phylogenetic analysis of all species of 
the genus Agonum that was based on numerous mor- 
phological characters recognised four subgeneric enti- 
ties: Agonum sensu stricto, Europhilus Chaudoir, 1859, 
Olisares Motschulsky, 1865 and Platynomicrus Casey, 
1920 (Liebherr and Schmidt 2004; Liebherr et al. 2005). 
Whereas some species are distinctive, others are difficult 
to identify as a consequence of intraspecific variations in 
setation, pronotal shape, body colouration and size and, 
therefore, typically require the examination of the geni- 
talia (Htrka 1996; Luff 2007). In the case of the Agonum 
species that are documented for Germany, the combina- 
tion of various characteristic traits, including, besides 
others, a well-developed median tooth at the mentum 
and the shape of the posterior angles of the pronotum, 
allows a correct identification (Schmidt, pers. commu- 
nication). In Europe, species related to Agonum viduum 
are most challenging in terms of identification, but this 
group has been carefully revised more than 25 years ago 
(Schmidt 1994). For Germany, 24 species have been re- 
corded so far (Trautner et al. 2014; Schmidt et al. 2016). 
As a consequence of many Agonum species being found 
in highly-endangered peat bogs or similar freshwater-as- 
sociated habitats, numerous species have become very 
rare or even threatened with extinction, for example Ago- 
num hypocrita (Apfelbeck, 1904) or Agonum munsteri 
(Hellén, 1935) (Schmidt et al. 2016). 

In this study, we present, as part of the on-going GBoL 
project, another step in generating a comprehensive 
DNA barcode library for the molecular identification of 
Central European ground beetle species, here focusing 
on the genus Agonum. The analysed sequence library 
included 23 species of Agonum and Oxypselaphus ob- 
scurus (Paykull, 1790) as outgroup. We generated 68 
new barcodes including some sequences of old pinned 
museum specimens and analysed a total number of 258 
DNA barcodes in detail. 


Material and methods 
Sampling of specimens 


Most of the sampled ground beetles (7 = 186, 78%) were 
collected between 1999 and 2015 using various sampling 
methods (i.e. hand collecting, pitfall traps). All beetles 
were stored in ethanol (96%). The analysed specimens 
were identified using the identification key provided in 
Schmidt (2006). It was also possible to generate DNA bar- 
codes from decade-old pinned ground beetles of the cara- 
bid collection of the Bavarian State Collection of Zoology, 
namely Agonum antennarium (Duftschmid, 1812) (n= 2) 
with an age of 87 years and A. impressum (Panzer, 1796) 
(n= 1, age: 52 years). In total, 68 new barcodes of 17 spe- 
cles were generated, including five species new to BOLD. 
For our analysis, we also included 175 DNA barcodes of 
four previous studies (Raupach et al. 2010: 17 specimens, 
5 species; Hendrich et al. 2015: 101 specimens, 13 spe- 
cies; Pentinsaari et al. 2014: 52 specimens, 12 species; 
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Figure 1. Representative images of analysed Agonum species. A. Agonum (Olisares) duftschmidi Schmidt, 1994; B. Agonum 
(Olisares) ericeti (Panzer, 1809); C. Agonum (Olisares) hypocrita (Apfelbeck, 1904); D. Agonum (Olisares) sexpunctatum (Linné, 
1758); E. Agonum (Agonum) marginatum (Linné, 1758); KF. Agonum (Agonum) muelleri (Herbst, 1784); G. Agonum (Europhilus) 
fuliginosum (Panzer, 1809); H. Agonum (Europhilus) piceum (Linné, 1758) and I. Agonum (Europhilus) thoreyi Dejean, 1828. Scale 
bars: 1 mm. All images were obtained from http://www.eurocarabidae.de. 
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Rulik et al. 2018: 5 specimens, 1 species), as well as 15 
public barcodes deposited in BOLD without correspond- 
ing publication (10 species). 

Most of the studied beetles were collected in Germany 
(n = 191, 73.9%), but for comparison, a number of spec- 
imens were included from Austria (7 = 10, 3.6%), Czech 
Republic (n = 1, 0.4%), Estonia (n = 1, 0.4%), Finland 
(n= 49, 18.9%), France (n= 2, 0.8%), Italy (n = 1, 0.4%), 
Montenegro (n = 2, 0.8%) and Sweden (n = 2, 0.8%). The 
number of analysed specimens per species ranged from 
one (A. gracilipes (Duftschmid, 1812), A. impressum 
(Panzer, 1796)) to a maximum of 31 for A. fuliginosum 
(Panzer, 1809). 


DNA barcode amplification, sequencing and 
data depository 


All laboratory operations were carried out, following 
standardised protocols for COI amplification and se- 
quencing (Ivanova et al. 2006; deWaard et al. 2008), 
at the Canadian Center for DNA Barcoding (CCDB), 
University of Guelph, the molecular labs of the 
Zoologisches Forschungsmuseum Alexander Koenig 
(ZFMK) in Bonn and the working group Systematics 
and Evolutionary Biology at the Carl von Ossietzky 
University Oldenburg, the latter two being located in 
Germany. Photos from each studied beetle were taken 
before molecular work started. One or two legs of one 
body side were removed for the subsequent DNA ex- 
traction which was performed using the QIAmp Tissue 
Kit (Qiagen GmbH, Hilden, Germany) or NucleoSpin 
Tissue Kit (Macherey-Nagel, Dtren, Germany), fol- 
lowing the extraction protocol. 

Detailed information about used primers, PCR ampli- 
fication and sequencing protocols can be found in a previ- 
ous publication (see Raupach et al. 2016). All purified PCR 
products were cycle-sequenced and sequenced in both 
directions by contract sequencing facilities (Macrogen, 
Seoul, South Korea or GATC, Konstanz, Germany), us- 
ing the same primers as used in PCR. Double-stranded 
sequences were assembled and checked for mitochondri- 
al pseudogenes (numts) by analysing the presence of stop 
codons and frameshifts, as well as double peaks tn chro- 
matograms with the Geneious Prime 2020.0.4 (https:// 
www.geneious.com) (Biomatters, Auckland, New Zea- 
land). For verification, BLAST searches (nBLAST, 
search set: others, programme selection: megablast) were 
conducted to confirm the identity of all new sequences 
based on already published sequences (high identity val- 
ues, very low E-values). 

Detailed voucher information, taxonomic classifications, 
photos, DNA barcode sequences, primer pairs used and 
trace files (including their quality) were uploaded to the pub- 
lic dataset “DS-BAAGO” (Dataset ID: dx.do1.org/10.5883/ 
DS-BAAGO) on the Barcode of Life Data Systems (BOLD; 
www.boldsystems.org) (Ratnasingham and Hebert 2007). 
Parallel to this, all new barcodes were deposited in Gen- 
Bank (accession numbers: MT520822—MT520889). 
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DNA Barcode analysis 


The complete dataset was analysed by using various 
approaches. First, the comprehensive analysis tools of 
the BOLD workbench were employed to calculate the 
nucleotide composition of the sequences and distribu- 
tions of Kimura-2-parameter distances (K2P; Kimura 
1980) within and between species (align sequences: 
BOLD aligner; ambiguous base/gap handling: pairwise 
deletion). In addition, all barcode sequences were sub- 
jected to the Barcode Index Number (BIN) analysis sys- 
tem implemented in BOLD that clusters DNA barcodes 
in order to produce operational taxonomic units that 
typically closely correspond to species (Ratnasingham 
and Hebert 2013). A threshold of 2.2% was applied for 
a rough differentiation between intraspecific and inter- 
specific distances, based on Ratnasingham and Hebert 
(2013). These BIN assignments on BOLD are constant- 
ly updated as new sequences are added, splitting and/or 
merging individual BINs in the light of new data. 

Second, maximum parsimony networks were con- 
structed with TCS 1.21, based on default settings 
(Clement et al. 2000) and implemented in the software 
package PopART v.1.7 (Leigh and Bryant 2015), in the 
case of species pairs with interspecific distances smaller 
than 2.2% and sharing identical BINs. Such networks al- 
low a better visualisation of the distances between close- 
ly-related species than classical tree topologies. 

Finally, all sequences were aligned using MUSCLE 
(Edgar 2004) and analysed using a neighbour-joining 
cluster analysis (NJ; Saitou and Nei 1987), based 
on K2P distances with MEGA 10.0.5 (Kumar et al. 
2018). Non-parametric bootstrap support values were 
obtained by re-sampling and analysing 1,000 pseudo- 
replicates (Felsenstein 1985). It should be noted that 
DNA barcodes do not aim to recover phylogenetic 
relationships (e.g. DeSalle and Goldstein 2019). Instead, 
the shown topology represents a graphical visualisation 
of DNA barcode distance divergences and putative 
species cluster. 


Results 


Overall, 252 DNA barcode sequences of 24 Agonum spe- 
cies were analysed, representing 96% of all document- 
ed species (n = 25) of this genus for Germany, except 
Agonum nigrum Dejean, 1828. A full list of the analysed 
species 1s presented in the supporting information (Suppl. 
material 1). Lengths of the analysed DNA barcode frag- 
ments ranged from 307 to 658 base pairs (bp). As is typi- 
cally known for arthropods, the DNA barcode region has 
a high AT-content (68%), with mean sequence compo- 
sitions for Adenosine (A) = 31%, Cytosine (C) = 16%, 
Guanine (G) = 16% and Tyrosine (T) = 37%. Intraspecif- 
ic K2P distances ranged from 0 to a maximum of 1.3% 
(Agonum ericeti (Panzer, 1809)), whereas interspecific 
distances had values between 0.16 and 7.07%. Lowest 
interspecific distances were found for Agonum micans 
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Figure 2. Maximum statistical parsimony networks of the species pairs. A. Agonum micans (Nicolai, 1822) (blue) and Agonum scitulum 
Dejean, 1828 (red); B. Agonum impressum (Panzer, 1796) (yellow) and Agonum sexpunctatum (Linne, 1758) (green) and C. Agonum dufts- 
chmidi Schmidt, 1994 (violet) and Agonum emarginatum (Gyllenhal, 1827) (light brown). Used parameters included default settings for 
connection steps, gaps were treated as fifth state. Each line represents a single mutational change, whereas small black dots indicate missing 
haplotypes. The numbers of analysed specimens (7) are listed, whereas the diameter of the circles is proportional to the number of haplo- 
types sampled (see given open circles with numbers). Scale bars: 1 mm. Beetle images were obtained from http://www.eurocarabidae.de. 
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(Nicolai, 1822) and Agonum scitulum Dejean, 1828) with 
a value of 0.16% and the same BIN (AAN9978). 

In total, unique BINs were revealed for 15 species (65%) 
and one BIN for two species for three species pairs (13%), 
namely Agonum micans (Nicolai, 1822)/Agonum scitulum 
Dejean, 1828), Agonum impressum (Panzer, 1796)/Agonum 
sexpunctatum (Linne, 1758) and Agonum duftschmidi 
Schmidt, 1994/Agonum emarginatum (Gyllenhal, 1827) 
(Fig. 2). Due to the fact that the numbers of unspecified nu- 
cleotides (“Ns”) exceeds more than 1% of their total length 
(Agonum antennarium (Duftschmid, 1812), 1 = 2) or only 
sequences < 500 bp were given (Agonum ericeti (Panzer, 
1809), n = 2), two species received no BIN assignment. 

The NJ analyses, based on K2P distances, revealed 
non-overlapping clusters for all analysed species (Fig. 3). 
In the case of species with more than two analysed spec- 
imens (n = 22), 82% (n = 18) of the studied species were 
characterised with bootstrap support values > 95%. A 
more detailed topology of all analysed specimens is pre- 
sented in the supporting information (Suppl. material 2). 


Discussion 


The results of this study highlight the efficiency of DNA 
barcodes for the determination of most German species 
of the genus Agonum, but also indicate a close relation- 
ship of a number of species with shared BINs, for ex- 
ample, Agonum scitulum (Dejean, 1828) and Agonum 
micans (Nicolai, 1822) (Suppl. material 1, Fig. 2). Hap- 
lotypes of both species are separated by a K2P distance 
of 0.16 and only three mutational steps (Fig. 2A). De- 
spite the fact that Agonum scitulum has been often over- 
looked and misidentified as Agonum micans as a result 
of flawed identification keys in the past (e.g. Paill 2010; 
Schmidt and Benedikt 2010; Brigi¢ et al. 2016), such a 
close relationship between these two species had, howev- 
er, not been expected before (see Liebherr and Schmidt 
2004). Agonum scitulum (Dejean, 1828) is a macropter- 
ous Western Palaearctic species with a discontinuous 
distribution from England to Romania, Croatia and the 
European part of Russia (Luff 2007; Paill 2010; Schmidt 
and Benedikt 2010; Brigi¢ et al. 2016). It is a very rare 
species that is typically found in marshes, fens and reed 
beds, syntopic with Agonum micans and Agonum fulig- 
inosum (Panzer, 1809) (Luff 2007; Paill 2010; Schmidt 
and Benedikt 2010). Thanks to a thorough review some 
years ago, an excellent key for the genus Agonum has 
been established and allows a reliable identification of 
both species, 1.c. the absence/presence of hairs on the 
dorsal site of the 3 tarsomere of the last walking leg 
(Schmidt 2006). Low molecular distances between both 
Species give evidence for an apparent recent separation, 
but the small number of studied specimens, as well as 
the only use of mitochondrial sequence data, does not 
allow more conclusions at the moment. Here, additional 
analysis combining a careful morphological analysis, as 
well as detailed fine-scaling nuclear sequence data, will 
provide more information. 
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In contrast to the previous species pair, a close rela- 
tionship between Agonum impressum (Panzer, 1796) and 
Agonum sexpunctatum (Linne, 1758), as well as Agonum 
dufischmidi Schmidt, 1994 and Agonum emarginatum 
(Gyllenhal, 1827), has been already suggested in the past 
(Liebherr and Schmidt 2004). The given barcode data 
support this hypothesis. Similar to the Agonum micans 
and Agonum scitulum, haplotypes of Agonum impressum 
(Panzer, 1796) and Agonum sexpunctatum (Linné, 1758) 
are separated by three mutational steps (K2P distance: 
0.36) (Fig. 2B), whereas the minimum K2P distance of 
Agonum dufischmidi Schmidt, 1994 and Agonum emargi- 
natum (Gyllenhal, 1827) has a value of 0.77, resulting in 
five additional mutational steps (Fig. 2C). Only the anal- 
ysis of additional specimens from different locations will 
show if these distances within both species pairs are sound. 

In comparison with other Central European ground 
beetles, for instance species of the genera Bembidion 
Latreille, 1802 (Raupach et al. 2016), Calathus Bonelli, 
1810 (Hendrich et al. 2015) or Notiophilus Duméril, 1806 
(Raupach et al. 2019), no interspecific distances > 2.2% 
or multiple BINs per species were found for the studied 
species and specimens. Highest intraspecific distances 
were identified for Agonum ericeti (Panzer, 1809) with a 
value of 1.3, Agonum muelleri (Herbst, 1784) (1.24) and 
Agonum thoreyi Dejean, 1828 (1.03) (Suppl. material 1). 
Nevertheless, all three species showed no conspicuous 
substructure for the analysed COI sequences yet. 

At this point it should be kept in mind that not all bee- 
tles were collected in Germany. For instance, all speci- 
mens of Agonum antennarium (Duftschmid, 1812) were 
sampled in Montenegro (see Suppl. material 2), whereas 
beetles of other species were sampled from different local- 
ities in Germany only (e.g. Agonum lugens (Duftschmid, 
1812) or Agonum micans (Nicolai, 1812)). As previous- 
ly mentioned, a number of species are very rare in the 
target area, but more abundant in adjacent countries. In 
most cases, the new DNA barcodes represent, however, 
the very first molecular data for these taxa and give im- 
portant impressions about their molecular diversity, even 
if they were not sampled in Germany. Nevertheless, it is 
also planned to analyse specimens for such species that 
were collected in Germany in the near future. 


Conclusion 


As central part of modern biodiversity research, DNA 
barcoding will become more and more prominent tn this 
research field. Therefore, comprehensive DNA barcode li- 
braries of important bio-indicators will become the back- 
bone of any applied analysis, for example, metabarcoding 
as well as eDNA studies. The present study demonstrates 
the successful identification of most Agonum species 
documented for Germany by using DNA barcodes. The 
dataset also shows a close relationship between various 
species and their putative recent origin. Only the analysis 
of additional species and specimens, however, will reveal 
if the observed patterns of distinct clusters persist. 
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Figure 3. Neighbour-joining (NJ) topology of the analysed ground beetle species, based on Kimura 2-parameter distances. Triangles 
show the relative number of individual’s sampled (height) and sequence divergence (width). Numbers next to nodes represent 
non-parametric bootstrap values > 75% (1,000 replicates). Beetle images were obtained from http://www.eurocarabidae.de except 
for Agonum monachum (Duftschmid, 1812) (Photo taken by Lars Hendrich). 
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Barcode analysis using the BOLD 
workbench 


Authors: Michael J. Raupach, Karsten Hannig, Jérome 
Moriniere, Lars Hendrich 

Data type: Data table 

Explanation note: Molecular distances are based on the 
Kimura 2-parameter model of the analysed specimens 
of the analysed species of the genus Agonum and 
Oxypselaphus obscurus (Paykull, 1790). Divergence 
values were calculated for all studied sequences, using 
the Nearest Neighbour Summary implemented in the 
Barcode Gap Analysis tool provided by the Barcode of 
Life Data System (BOLD). Align sequencing option: 
BOLD aligner (amino acid-based HMM), ambiguous 
base/gap handling: pairwise deletion. ISD = intraspe- 
cific distance. BINs are based on the barcode analysis 
from 16-04-2020. Species pairs with interspecific dis- 
tances < 2.2% are marked in bold. As a consequence of 
that the number of unspecified nucleotides (“Ns”) ex- 
ceeds more than 1% of the total length of the sequence, 
barcodes of Agonum antennarium (Duftschmid, 1812) 
and A. ericeti (Panzer, 1809) received no BIN. Country 
codes: D = Germany, A = Austria, B = Belgium, BG = 
Bulgaria, CZ = Czech Republic, EST = Estonia, FIN 
= Finland, F = France, I = Italy, RO = Romania, SK 
= Slovakia, SLO = Slovenia, and SW = Switzerland. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0). The Open Database License 
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(ODbL) is a license agreement intended to allow us- 

ers to freely share, modify, and use this Dataset while 

maintaining this same freedom for others, provided 

that the original source and author(s) are credited. 
Link: https://do1.org/10.3897/dez.67.56163.suppl1 


Supplementary material 2 
Neighbour-joining topology 


Authors: Michael J. Raupach, Karsten Hannig, Jérome 
Moriniere, Lars Hendrich 

Data type: Neighbour-joining topology 

Explanation note: Neighbour-joining topology of all an- 
alysed carabid beetles based on Kimura 2-parameter 
distances. Specimens are classified using ID numbers 
from BOLD and species name. Numbers next to nodes 
represent non-parametric bootstrap values (1,000 rep- 
licates, in %). Three beetles of a previous publication 
that were identified as Agonum ericeti (Panzer, 1809) 
(see Pentinsaari et al. 2014) were incorrectly deter- 
mined. A careful re-inspection revealed these speci- 
mens as Agonum sexpunctatum (Linné, 1758). 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://doi.org/10.3897/dez.67.56163.suppl2 
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